Robust Patterns in Food Web Structure 
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We analyze the properties of seven community food webs from a variety of environments — 
including freshwater, marine-freshwater interfaces and terrestrial environments. We uncover quan- 
titative unifying patterns that describe the properties of the diverse trophic webs considered and 
suggest that statistical physics concepts such as scaling and universality may be useful in the de- 
scription of ecosystems. Specifically, we find that several quantities characterizing these diverse food 
webs obey functional forms that are universal across the different environments considered. The 
empirical results are in remarkable agreement with the analytical solution of a recently proposed 
model for food webs. 
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In natural ecosystems species are connected through 
trophic relationships [|], |[ ^ defining intricate networks 
l@' il H' 0' so-called food webs. Understanding the 
structure and mechanisms underlying the formation of 
these webs is of great importance in ecology ||] . For this 
reason, much research has been done in constructing em- 
pirical webs and uncovering unifying patterns describing 
their structure ||, ^. However, in the last decade the con- 
struction of larger and more complete food webs clearly 
indicated that the previously reported unifying patterns 
do not hold for the new webs @ |ll|. Indeed, the com- 
plexity of the new webs has rendered quite difficult the 
challenge to obtain quantitative patterns that substitute 
the old ones. 

Here, we analyze the properties of seven detailed com- 
munity food webs from a variety of environments — 
including freshwater habitats, marine-freshwater inter- 
faces and terrestrial environments. Remarkably, we un- 
cover quantitative unifying patterns that describe the 
properties of the diverse trophic webs considered and cap- 
ture the random and non-random aspects of their struc- 
ture. Specifically, we find that several quantities — such 
as the distributions of number of prey, number of preda- 
tors, and number of trophic links — characterizing these 
diverse food webs obey robust functional forms that de- 
pend on a single parameter, the linkage density z. 

In our analysis, we use results obtained for complex 
networks and for a recent model of food web forma- 
tion, the "niche model" of Ref. We first describe 
the theoretical model and its predictions: An ecosystem 
with S species and L trophic interactions between these 
species, defines a network with S nodes and L directed 
links. In the niche model, one first randomly assigns 
species i ^ 1 , . . . , to "trophic niches" which are 
mapped into the interval [0,1]. A species i is character- 
ized by its niche parameter and by its list of prey. Prey 
are chosen according to the following procedure: species 
i preys on the species j with niche parameters rij inside a 
segment of length xrii centered in a position chosen ran- 
domly inside the interval [xni/2,ni]. Here, < a; < 1 
is a random variable with probability density function 
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H. The values of the parameters 



Pxix) = b{l _ 

b = {S^/2L — 1) and S determine the linkage density 
z = L/S of the food web, and the directed connectance 
L/S'^, which is a measure of the fraction of the actual 
number of trophic links as compared to the maximum 
possible number 

In the limit of large web sizes {S ^ 1) and small 
connectances [L/S"^ <C 1), one can derive analytical 
expressions for the distribution of number of prey k 
|12[ . We consider the cumulative distribution Ppi.oy(A:) = 
X]fe'>fe-Ppicy(fc') because it is less noisy than the proba- 
bility function Pproy(^)- We obtain 
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where Ei(x) is the so-called exponential- integral func- 
tion Equation (|l|) predicts that the distribution of 
number of prey decays exponentially for large k. 

Also in the limit of large web sizes and small con- 
nectances, one can derive analytical expressions for the 
distribution of number of predators m 12 1. We obtain. 
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where j{m + l,z) is the so-called "incomplete gamma 
function" fl^ . To gain intuition about the functional 
form (|^), note that Ppredim) is approximately a step func- 
tion: It is constant for m < z, and then it decays with 
a Gaussian tail for m « 2z |jl^. It follows then that 
the cumulative distribution Pprcdi'Ti) decreases linearly 
as 1 — m/ z for m < z and decays as the error function 
||l3| for m w 2z. 

Next, we analyze the empirical data for seven food 
webs with 25 to 92 trophic species. These webs have 
linkage densities 2.2 < z < 10.8, and connectances in 
the interval 0.06-0.31 |]l|. We first investigate the dis- 
tributions of number of prey and number of predators. 
Figures |l|a,b compare the cumulative distributions of the 
number of prey and number of predators for species in 
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FIG. 1: Cumulative distribution (a) Pproy of number of prey 
k, and (b) Pprod of number of predators m for the St. Martin 
Island web. The data agrees well with the analytical predic- 
tions of Eqs. (0)-(^) — indicated by the solid lines, without 
any free parameters for fitting as z is determined empirically. 



the St. Martin Island web with our analytical pre- 
dictions, and suggest that these distributions are well ap- 



proximated by Eqs. (|1| 
for fitting. Equations (h) 




) without any free parameters 
and the results of Fig. |l| 



suggest the possibility that -Pprcy and Pprcd obey univer 



sal functional forms that depend only on z. 

Indeed, Eq. (|^) predicts that Ppfey(fc) depends only on 
k/2z. So, we plot in Figs. |^a,c the cumulative distribu- 
tions Ppi-ey(fc) versus the scaled variable k/2z for the food 
webs and find that the data collapse onto a single curve, 
supporting the possibility that Pproy obeys a universal 
functional form |l5|| . 

The scaling of Ppred("T') is not as straightforward. 
Equation (H) indicates that "true" scaling holds only for 
m/2z < 1/2, while for larger values of m/2z there is a 
Gaussian decay of the probability function with an ex- 
plicit dependence on z. However, the decay for m > 2z is 
quite fast and, to first approximation, not very relevant. 
Thus, we plot Pprcd (jTi) versus the scaled variable to/2z 
for the food webs and indeed find a collapse of the data 
onto a single curve for m/2z < 0.7 (Figs. ||b,d) [l^ . 

Figure || supports the strong new hypothesis that the 
distributions of number of prey and number of predators 
follow universal functional forms. To improve statistics 
and better determine the specific functional form of these 
distributions, one may pool the scaled variables, k/2z 
and m/2z, from the different webs |l6| into single 
distributions, Pprey and Pprod, respectively. Figures ||a,b 
show the cumulative distributions of scaled number of 
prey and scaled number of predators. Note that the dis- 
tributions are well approximated by Eqs. ([^)-(|^) even 
though there are no free parameters to fit in the analyt- 
ical curves. These results are analogous to the finding of 
scaling and universality in physical, chemical and social 
systems Q. 

Figure He plots the probability densities for the distri- 
bution of number of prey and number of predators. It 
is visually apparent that both distributions are difi^erent. 
This is confirmed by the Kolmogorov-Smirnov test which 
rejects the null hypothesis at the p < 0.001 level. The 
distribution of number of prey decays exponentially, and 
the distribution of number of predators is essentially a 
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FIG. 2: We test the "scahng hypothesis" Q that the dis- 
tributions of number of prey (predators) have the same func- 
tional form for different food webs, (a) Cumulative distri- 
bution Pproy of the scaled number of prey k/2z for all the 
webs except Ythan |l6j . The solid line is the prediction of 
Eq. (|l|). The data "collapses" onto a single curve that agrees 
well with the analytical results, (b) Cumulative distribution 
Pprod of the scaled number of predators m/2z for all the webs 
but Ythan . The solid lines are the analytical predic- 

tions of Eq. (El) for the extremal values of z in the empirical 
data. Semi-logarithmic plot of the scaled distributions of (c) 
number of prey, and (d) number of predators. 



step function with a fast decay. 

One can perform a similar analysis for the distribution 
Piink of the number of trophic links r = k + m. As for 
number of prey or number of predators, the data from 
the different webs, upon the scaling r/2z, collapse onto 
a single curve, further supporting the hypothesis that 
scaling holds for food web structure. To better deter- 
mine the specific functional form of piink(?'), we pool the 
scaled variables, r/2z, from all webs except Ythan into 
a single distribution (Fig. ^). We find that punkir) has 
an exponential decay for r/2z ^ 1, in agreement with 
our theoretical calculations. Therefore, there is a char- 
acteristic scale for the linkage density, i.e. food webs do 
not have a scale-free structure, in contrast to reports in 
recent studies of food- web fragility |l9|| . 

Next, we test if the scaling hypothesis suggested by 
the analysis of distribution of trophic links also applies 
to other quantities characterizing food web structure. We 
consider two quantities with ecologic implications: (i) the 
average trophic distance d between species Q (which is 
the number of species needed to trophically connect two 
given species), and (ii) the clustering coefficient C (which 
counts the fraction of species' triplets that form fully- 
connected triangles). The latter relates to the compart- 
mentalization in an ecosystem while the former relates to 
typical food-chain length. 
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ber of prey fc/22 for the pooled webs (all except Ythan). The 
solid line is the analytical prediction (^) , and the dashed line 
is a numerical simulation of the niche model fl | with S — 244 
(the size of the pooled data) and z — 7.5 (the average de- 
gree for the pooled webs), (b) Cumulative distribution Pprcd 
of the scaled number of predators m/2z for the pooled webs. 
The solid line is the analytical prediction (^) for the case 
z = 7.5, and the dashed line is a numerical simulation of the 
niche model |Ql with S = 244 and z = 7.5. (c) Comparison 
of the probability density functions of the scaled number of 
prey and number of predators. It is visually apparent that the 
two distributions have distinct functional forms, (d) Proba- 
bility density function of the number of trophic interactions 
per species r = k + m pooled for all webs except Ythan. The 
solid line is obtained by numerically convolving the distribu- 
tions (0)-(^ while the dashed line is obtained by numerical 
simulations of the niche model in the limit of large web sizes 
and small connectances (we use S = 1000 and z = 5) — the 
same limit for which the analytical curves where derived |^2| . 
The tail of the distribution decays exponentially, indicating 
that food webs do not have a scale-free structure. 
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FIG. 4: (a) Scaled average trophic distance d between species 
versus linkage density z. We compare the data with the 
numerical simulations of the niche model for web sizes 
S = 100,500,1000 (thin solid lines). We find a logarith- 
mic increase of the average distance with web size S, in good 
agreement with the model predictions. We also compare our 
results with the prediction for a random graph with the same 
linkage density as the webs studied (dashed line). The log- 
arithmic dependence of d{S) agrees with the expectation for 
a random graph, however, the coefficients of the logarithmic 
increase differ from the predicted values indicating that food 
webs have a more complex structure than that predicted by a 
random graph, (b) Double-logarithmic plot of the clustering 
coefficient C versus the scaled web size S/2z. We compare 
the data with numerical results for the niche model [Q| for 
three values of the linkage density in the empirically-relevant 
range {z = 2.5, 5, 10). We find that the clustering coefficient 
of the food webs is inversely proportional to the web size S, 
in good agreement with the model predictions and with the 
asymptotic behavior predicted for a random graph M. 



In Fig. |ja, we compare our numerical results for the 
average trophic distance d for the niche model with 
the values calculated for the food webs analyzed. We find 
that d increases with web size as log S both for the model 
and for the data. This logarithmic increase is the ex- 
pected behavior for a random graph; however, the slopes 
measured for the data and the model are different from 
the value predicted for a random graph 0, suggesting 
that there is a degree of "order" to the connectivity of 
the food web which may encode the mechanisms of food 
web assembly. Remarkably, this characteristic of the em- 
pirical food webs appears to be captured by the niche 
model The results of Fig. ^ also support the scal- 
ing hypothesis and suggest that the average distance in 
a food web may also follow a unique functional form for 
different food webs. 

Figure 0b shows our results for the clustering C of 



the food webs studied and for the niche model |^]. We 
find that the data is well approximated by the model 
predictions, and that C decreases to zero as l/S* as web 
size S increases. 

The major finding of this study is the uncovering of 
unifying quantitative patterns characterizing the struc- 
ture of food webs from diverse environments. Specifi- 
cally, we find that the distributions of the number of prey, 
number of predators and number of links of most of the 
best studied food webs seem to collapse onto the same 
curves after rescaling the number of links by its average 
number z. Remarkably, the corresponding curves are in 
agreement with the analytical predictions of the niche 
model. Therefore, these distributions can be theoreti- 
cally predicted merely by knowing the food web's link- 
age density z, a parameter readily accessible empirically. 
Regularities such as these are interesting as descriptors 
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of trophic interactions inside communities because they 
may enable us to make predictions in the absence of high- 
quahty data, and provide insight into how communities 
function and are assembled. 

Our results are of interest for a number of other rea- 
sons. First, food webs do not have a scale free distri- 
bution of number of links (total, incoming or outgoing) . 
This is surprising since one could expect most species to 
try to prey on the most abundant species in the ecosys- 
tem (an "abundant-get-eaten" kind of mechanism) . Such 
a preferential attachment would lead to a scale-free dis- 
tribution of links; instead, we find a single-scale distri- 
bution, suggesting that species specialize and prey on a 
small set of other species. Second, the results of Figs. ||a,b 
support the scaling hypothesis and indicate that there 
is very little, if any, compartmentalization in ecosystems 
po| , suggesting the possibility that ecosystems are highly 
interconnected and that the removal of any species may 
induce large disturbances. Third, the structure of food 
webs is different from many other biological networks in 
two important aspects: the links are uni-directional and 
the in- and out-degree distribution are different. These 
two facts are a result of the directed character of the 
trophic interactions and of the asymmetry it creates. In- 
terestingly, the niche model captures this asymmetry in 
its rules, which may explain its success in explaining the 



empirical results. 

Our findings are surprising for two reasons: (i) they 
hold for the most complete food webs studied, in contrast 
to previously reported patterns and (ii) they sup- 
port the possibility that fundamental concepts of mod- 
ern statistical physics such as scaling and universality 
JlSf — which were developed for the study of inanimate 
systems — may also be applied in the study of food webs — 
which comprise animate beings. Indeed, our results are 
consistent with the underlying hypothesis of scaling the- 
ory p^ , i.e., food webs display universal patterns in the 
way trophic relations are established despite apparently 
"fundamental" differences in factors such as the environ- 
ment (e.g. marine versus terrestrial), ecosystem assem- 
bly, and past history. 
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